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ABSTRACT 

We provide a data set of emission and recombination coefficients of hydrogen using a re- 
distribution of electron energies rather than the more traditional Maxwell-Boltzmann 
(MB) distribution. The data are mainly relevant to thin and relatively cold plasma 
found in planetary nebulae and H II regions. The data set extends the previous data 
sets provided by Storey and Hummer which were computed using a MB distribution. 
The data set, which is placed in the public domain, is structured as a function of 
electron number density, temperature and k . Interactive fortran 77 and C-| — b data 
servers are also provided as an accessory to probe the data and obtain Lagrange- 
interpolated values for any choice of all three variables between the explicitly computed 
values. 

Key words: atomic data - atomic processes - radiation mechanisms: general - radi- 
ation mechanisms: non-thermal - ISM: abundances - planetary nebulae: general. 


1 INTRODUCTION 

There are many studies related to the recombination of hy- 
drogen and hydrogenic systems, the most comprehensive 


being those of Hummer & Storey (1987); Storey & Hum- 


11988) and Storey & Hummer (1995 b However, all 


the past work was based on a Maxwell-Boltzmann distri- 
bution of electron energies. There has been, and still is a 
general consensus that this is the appropriate distribution 
for thin nebular plasmas. However, this was disputed in the 


past (Hagihara 1944) where considerable deviations from the 


thermodynamic equilibrium on which the MB relies were 
claimed although this claim was later discounted |Bohm & 
Aller| 1947 1. 


There has been a recent revival (Nicholls et al 20121 


of the proposal of a non-thermal electron energy distribu- 
tion in planetary nebulae and H II regions in the light of 
the long standing problem in nebular physics of the con- 
tradiction between the results for elemental abundance and 
electron temperature as obtained from the optical recombi- 
nation lines (ORL) and those obtained from the collision- 
ally excited lines (CEL). According to the recent proposal, 
the ORL-CEL discrepancy problem can be resolved by as- 
suming a non-MB electron distribution, specifically the k- 
distribution. There have been a few recent attempts to as- 
sess the merit of this suggestion ( Sochi|2012 Storey & Sochi 


20131 [20l4] |Zhang et at^OU I using spectroscopic means to 


directly sample the free electron energy distribution. With 
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the exception of Storey & Sochi (2014) they are all inconclu- 


sive, in the sense that the data do not differentiate between 
a single K-distribution and other models, such as one with 
two MB components at different temperatures. |Storey fc] 


Sochi (20141 do, however argue that, with a high degree of 


certainty, the Balmer line and continuum spectrum of the 
extreme planetary nebula Hf 2-2 cannot be modeled with 
a single /-e-distributed electron energy distribution, while it 
can be with a model comprising two MB distributions. It 


should be noted however that Zhang et al (2014) analyze 


the same spectra of the same object and conclude that ei- 
ther model can adequately model the spectrum. We return 
to this apparent contradiction below. 

Both |Storey fe Sochi| ( |2014| ) and |Zhang et al\ ( |2014| ) 
model the Balmer line and continuum spectrum with MB 
and k distributions. Typically the continuum intensity is 
modeled relative to one of the high Balmer lines, chosen to 
be close in wavelength to the Balmer edge and apparently 
unblended. The continuum spectrum is relatively easy to 
model with an arbitrary electron energy distribution but the 
modeling of the high Balmer line intensities requires a full 
treatment of the collisional-radiative recombination process 
as a function of k as well as the usual density and tempera- 


ture variables. Storey & Sochi (2014) made such a calcula- 


tion in Case B of Baker & Menzel ( 1938 1 and presented some 


results for the line which they used to normalize intensities, 
Hll. Here we publish the full results from those calculations. 


We note that Zhang et al (2014) also use Hll for normal- 
ization but they rely on an approximate treatment taken 


from Nicholls et al (2012) which applies a K-dependent scal- 
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ing function to the Hll emission coefficient calculated with 
a MB distribution. Storey & Sochi (20141 show that this 
approximation is poor for the very low values of k that are 
required to model the Hf 2-2 spectrum, which may explain 


why Zhang et al (2014) reach a different conclusion to Storey 


& Sochi (2014) in the case of Hf 2-2. 


The results are provided in two text files where the emis- 
sion and recombination coefficients are given as a function 
of electron number density N e , electron temperature T e and 
k. We also provide interactive data servers, in the form of 
fortran 77 and C++ codes, for mining the data and obtain- 
ing interpolated values in the three variables between the 
explicitly computed values. In section [2] we give a brief the- 
oretical background about the atomic computational model 
used to generate the data, while in section [3] we explain the 
structure of the data set and provide general instructions 
and clarifications about how it should be probed and used. 
Section [4] contains general conclusions and discussions. 


2 ATOMIC MODEL 


The fundamental quantity required for the calculation of 
the hydrogen recombination line intensities is the coefficient 
for recombination to a state nl of H. For a free electron 


energy distribution f(E), this is given by (e.g. Storey & 
Sodul png) 


y(nl) = 


Ri 


j(7ll) 
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where R is the Rydberg energy constant, lo + and ai(nZ) are 
the statistical weights of the initial and final states respec- 
tively, v is the frequency of the emitted photon, E is the 
energy of the free electron, a(u, nl) is the cross-section for 
photoionization which is the inverse process to recombina- 
tion, and the other symbols have their standard meanings. 
The calculation of the recombination line emission coeffi- 
cients in a full collisional-radiative treatment has been de- 
scribed by |Hummer fe Storey| ( fT~987[ ) and |Storey fc Hummer| 
( 1995 ) and we use the same methods here. 


Traditionally, f(E) is the Maxwell-Boltzmann distribu- 
tion function. Here we use instead the k distribution which 
is given by (Vasyliunas 1968, Summers fe Thorne|199 1 1 


f K (E,T K ) = 

2 Ve r(/c + 1) ^ __e__Y ( k+ 

V7r(fcT K )3 ( K _ f)lr(K_ |) V + (K- l)kTj 

(2) 

where c is a parameter defining the distribution, Y is the 
gamma function of the given arguments, and T K is a temper- 
ature characteristic to the particular distribution. For suffi- 
ciently large k, f K (E,T K ) tends to the Maxwell-Boltzmann 
distribution. 

For energetically low-lying states of H the dominant 
processes are recombination and radiative decay. For higher 
states, collisional processes become important, with Z- 
changing collisions being the most frequent. In our calcu- 
lations the nl states of hydrogen are assumed degenerate 


with respect to l and in this case the dominant processes 
that change l are collisions not with electrons but with H + , 
He + and He ++ ions. In the results described here we retain a 
Maxwell-Boltzmann distribution for these heavier particles. 
We also retain a Maxwell-Boltzmann distribution for the 
processes that change energy and n, which are dominated by 
collisions with electrons. Consequently the emissivities that 
we calculate should be treated with caution for the high-n 
states for which Z- and n-changing collisional processes be- 
come important. The boundary of this region is primarily 
a function of the electron density, being at approximately 
n = 100,75,50,30 and 20 for densities of 10 2 , 10 3 , 10 4 , 10 5 
and 10 6 cm -3 respectively. 

The rate coefficients for Z-changing collisions used to ob- 
tain the above boundary n values were calculated using the 
theory described by |Pengelly fe Seaton| (|1964[) . Vrinceanu 


and co-workers have published a series of papers ( Vrinceanu 


2005 Vrinceanu et al 2012 2014) on electron and pro- 


ton induced collisions with Rydberg states of hydrogen. In 


Vrinceanu et al (2012) they state that the rate coefficients 


for proton induced AZ = 1 transitions are overestimated by 
the theory of |Pengelly &; Seaton| ( |1964[ ) by about an order of 
magnitude. The calculation of the rate coefficient depends 
upon an integration of the probability for an Z-changing 
transition over the impact parameter of the incident par- 
ticle, assumed to travel on a straight-line trajectory. It is 
well known that this integral is divergent for AZ = 1 transi- 
tions in a quantum mechanical treatment. The approximate 
treatment of the transition probability by Pengclly & Seaton] 
(1964) converges on the quantum mechanical result at large 
impact parameter as illustrated in Figure 1 of | Vrinceanu eZ| 
[aZ (2012 1. Pengelly & Seaton (19641 introduce a cut-off at 
large impact parameter to remove the divergence based on 
collective effects in the plasma or the finite lifetime of the 
Rydberg state. The semi-classical approach of | Vrinceanu et\ 
[aZ (20121 does not correctly replicate the quantum behav- 
ior at large impact parameter with the probability instead 
falling discontinuously to zero at a finite and relatively small 
value of the impact parameter. The missing contribution 
from large impact parameter is the origin of the order of 
magnitude difference they report between their results and 
those of Pengelly & Seaton ( 1964). We see no reason to pre- 
fer their semi-classical result over the quantum treatment at 
large impact parameters and therefore consider the|Pengelly| 


& Seaton ( 1964 1 results to be more reliable. 


3 DATA 

The K-dependent emission coefficients, e(ZV e ,T e , k) are pro- 
vided in a single text file called ‘elbk.d’. The energy emit- 
ted per unit volume per unit time is then N e N + e(N e , T e , k) 
where ZV+ is the H + number density and where all quantities 
are in cgs units. The structure of this file is explained in the 
following bullet points: 

• The first row of this file contains (in the following order) 
the number of N e values, the number of T e values and the 
number of k, values for which data are provided. 

• The 9 values of ZV e are specified by logioA e = 2. 0(0. 5)6.0. 

• The 16 values of T e are specified by logioTe = 2. 0(0. 2)3. 8, 
3.9(0. 1)4.4. 
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The 44 values of re are specified by logiore = 0.20(0.01)0.30, 
0.35(0.05)1.0, 1.1(0. 1)2.0, 2.2(0.2)3.0, 3.5, 4.0, 5.0, 6.0. 

The data therefore consist of 6336 blocks (= 9 x 16 x 44). 
The first row of each block contains information about the 
block which consists of the following: 

Z loglofVe logloTe logiore B Tlmin Tbmax 
where Z = 1 is the atomic number of hydrogen, ‘B’ refers to 
Case B and rimin and Umax are the minimum and max- 
imum upper state principal quantum numbers for which 
emission coefficients are tabulated. Each block therefore con- 
tains 4850 (= |n m (n m — 1) — 1) entries. These 4850 entries 
are arranged in 607 rows and hence the total number of rows 
in each block is 608. 

Of the three variables, N e T e and re, the fastest varying is N e 
followed by T e followed by re, and hence the ordinal number 
of a block is given by: 

Ob = On + (O t - 1)9 + (O k - 1)144 (3) 

where Ob, On, Ot and O k are the ordinal numbers of block, 
N e value, T e value and re value respectively. For example 
the ordinal number of the block for logioAL = 4 (On = 5), 
logioT e = 2.6 (Ot = 4) and logiore = 0.27 (O k = 8) is: 

O b = 5 + (4 - 1)9 + (8 - 1)144 = 1040 

and hence it starts on row 631714 (= 608(Os — 1) + 2) and 

ends on row 632321 (= 6080s + 1). 

The 4850 values of emission coefficients in each block are 
arranged for transitions from upper levels n u to lower levels 
ni with n u in descending order from n m , and rq in ascending 
order from 1 to (n u — 1), and hence the ordinal number for 
a transition tr(n u , n{) is given by 

O tr = ni + ^(n m - n u )(nm + n u ~ 1) (4) 

Emission coefficients to n = 1 are not calculated for Case B 
and hence are set to zero. 

A second file named ‘tlbk.d’ contains the hydrogen to- 
tal recombination coefficients in Case B, a(N e ,T e , k), and 
the total recombination coefficients to the 2s state of hydro- 
gen, ct2s(N e ,T e , re ), such that the number of recombinations 
per unit volume per unit time is N e N+a in cgs units. This 
file contains 12672 entries arranged in 1584 rows. The first 
half (6336) of these entries are the total recombination co- 
efficients of hydrogen while the second half are the total re- 
combination coefficients to the 2s state. The entries in each 
one of these two blocks correspond to the 6336 values of 
physical conditions (i.e. various combinations of N e , T e and 
re) positioned according to equation [3] 


provide easy access to the data with Lagrange-interpolated 
values in all three variables. 
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4 CONCLUSIONS AND DISCUSSIONS 


In this paper we computed atomic emission and recombi- 
nation coefficient data for hydrogen with electron energies 
described by a re-distribution. The atomic model used in the 
computation of these data uses the techniques described by 


Hummer & Storey (1987) and Storey & Hummer (19951. 


The data, which are placed in the public domain, span 
ranges of electron number density, temperature and re useful 
for modeling and analyzing plasmas such as those found in 
planetary nebulae and H II regions. Interactive data servers 
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